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Abstract. We consider the nonlinear Schrodinger equation on a compact manifold near an elliptic peri- 
odic geodesic. Using a geometric optics construction, we construct quasimodes to a nonlinear stationary 
problem which arc highly localized near the periodic geodesic. We show the nonlinear Schrodinger evolu- 
tion of such a quasimode remains localized near the geodesic, at least for short times. 



1. Introduction 

In this note, we study the Nonlinear Schrodinger Equation (NLS) on a compact, Riemannian manifold 
with a periodic elliptic (or stable) geodesic, which we define and discuss in more detail in Appendix |^ 
Specifically, we study solutions to 

where is the standard (negative semidcfinitc) Laplace-Beltrami operator and the solution is of the form 

To solve the resulting nonlinear elliptic equation, which can be analyzed using constrained variations, we 
will use Fermi coordinates to construct a nonlinear quasimode similar to the one presented in jThoOSb] for 
an arbitrary manifold with a periodic, elliptic orbit. To do so, we must analyze the metric geometry in a 
neighbourhood of a periodic orbit, for which we use the presentation in |Gra90j . For further references to 
construction of quasimodes along elliptic geodesies, see the seminal works Ralston |Ral82j . Babic |Bab68j . 
Guillemin-Weinstein |G W76| . Cardoso-Popov |CP02j . as well as the thorough survey on spectral theory 
by Zelditch jZelOSj . For experimental evidence of the existence of such Gaussian beam type solutions in 
nonlinear optics as solutions to nonlinear Maxwell's equations, see the recent paper by Schultheiss et al 
[SBS+IO] , 



Henceforward, we assume that (A/, g) is a compact Riemannian manifold of dimension d > 2 without 
boundary, and that it admits an elliptic periodic geodesic F C M . For /? > 0, we introduce the notation 

Ufi^{xe M : distg(a;,r) < P}. 

Let L > be the period of F. Throughout the paper, each time we refer to the small parameter /i > 0, 
this means that h takes the form 

(1.2) h=^, 

for some large integer TV e N. 

Our first result is the existence, in the case of a smooth nonlinearity, of quasimodes which arc highly 
concentrated near the elliptic periodic geodesic F. More precisely : 
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Theorem 1.1. Let {M,g) be a compact Riemannian manifold of dimension d>2, without boundary and 
which has an elliptic periodic geodesic. Let p be an even integer, let s > and assume that 

p{^~s)<l. 

Then for ft- <C 1 sufficiently small, and any S > 0, there exists iph{x) G H°°{M) satisfying 

(i) Frequency localization : For all r>Q, there exist C'i,C2 independent of h such that 

W-'' <\\'Ph\\H'-iM)<C2h'''\ 

(ii) Spatial localization near T : 

\\fh\\H-{M\U^i^2-s) 0{h°°)\\iph\\H^M)- 

(iii) Nonlinear quasimode : There exists X{h) e R so that iph satisfies the equation 

The quasimode iph and the nonhnear eigenvalue X{h) wih be constructed thanks to a WKB method, 
this win give a precise description of these objects (see Section \7.'3\\ . For instance, X{h) reads X{h) = 
- Eoh-^ + o(/i-i), for some Eq € R. 

As a consequence of Theorem ll.il and under the same assumptions, we can state : 

Theorem 1.2. Consider the function iph given by Theorem \Ll[ There exists cq > 0, such that if we denote 

by Th = CQh^^~^ ln( — ), the solution of the Cauchy problem (jl.ip with initial condition (fh satisfies 

\\uh\\L--ao,nhH-{M\U^^^2-5)) ^ Ch'^'^+^^^^. 

To construct a quasimode with 0{h°°) error, we need that the nonhnearity is of polynomial type, that 
is p € 2N. However, in the special case where d = 2 and s = 0, we are able to obtain a result for any 
< p < 4, which is a weaker version of Theorem 11.21 



Theorem 1.3. Let {M,g) be a compact Riemannian surface without boundary which admits an elliptic 
periodic geodesic and let < p < 4. Then for /i <C 1 sufficiently small, and any 5,e > 0, there exist 
(phix) € H°°{M) and v > Q satisfying 

(i) Frequency localization : For all r > 0, there exist Ci,C2 independent of h such that 

Clh-- < yh\\HnM)<C2h-''. 

(ii) Spatial localization near T : 

\\y^h\\L^M\U^i,2-5) = 0{h°°)\\iph\\L2(M), 

such that the corresponding solution Uh to (|l.ip satisfies 

for T/i = HP in the case p £ (0, 4)\{1}, and T^ ~ h^^^ in the case p = 1. 

Remark 1.1. In our current notations, the critical exponent for is 

_ 4 _ 4 
d-2 ^ n- 1' 

where n is the dimension of the geodesic normal hypersurface to T. The L^ critical exponent is p = 4/d, 
and the H^/'^ critical exponent is p ^ 4/((i — 1). 
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In general we construct highly concentrated solutions along an elliptic orbit, which is effectively a d — I 
dimensional soliton. Since stable soliton solutions exist on K'' for monomial Schrddinger of the type in 
()l.ip precisely for P < 2' numerology matches well that required for the expected local existence of 
stable nonlinear states on these lower dimensional manifolds. 

In addition, solving on compact manifolds with the power p = A/{d— 2) plays an important role 

in the celebrated Yamabe problem, see jSY94| . The authors hope that the techniques here can give insight 
into the energy minimizers of such a problem related to the geometry. 

Remark 1.2. In [CMlOj . the existence of nonlinear bound states on hyperbolic space, M'^ , was explored. 
In that case, it was found that the geometry made compactness arguments at oo rather simple, hence the 
only driving force for the existence of a nonlinear bound state was the local behavior of the nonlinearity. 
Work in progress with the second author, third author, Michael Taylor and Jason Metcalfe is attempting 
to generalize this observation to any rank one symmetric space. In this note, we find nonlinear quasimodes 
based purely on the local geometry, where here the nonlinearity becomes a lower order correction, which is 
a nicely symmetric result. 

Remark 1.3. In a private conversation with Nicolas Burq, he has pointed out that in settings where one 
assumes radial symmetry of the manifold, it is possible to construct exact nonlinear bound states. 

Notations. In this paper c, C denote constants the value of which may change from line to line. We use 
the notations a ^ b, a <b if ^b < a < Cb, a < Cb respectively. 

Acknowledgments. P. A. was supported by an NSF joint institutes postdoctoral fellowship (NSF 
grant DMS-0635607002) and a postdoctoral fellowship of the Foundation Sciences Mathcmatiques de Paris. 
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2. Outline of Proof 

The proof of Theorem 11.21 is to solve approximately the associated stationary equation. That is, by 
separating variables in the t direction, we write 

ilj{x,t) = e-'^*u{x), 

from which we get the stationary equation 

{X + Ag)u^a\u\Pu. 

For /i ^ 1 of the form (|1.2p , the construction in the proof finds a family of functions 

Uhix) = h-^^'-^^/^gih-^/^x) 

such that g is rapidly decaying away from F, C°°, g is normalized in L^, and 

(A + Ag)uh = a-\uhfuh + Q{uh), 

where A ~ and where the error Q{uh) is expressed by a truncation of an asymptotic series similar to 
that in |Tho08b| and is of lower order in h. 
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The point is that Theorem 1 1.2 1 (as well as Theorem 1 1.3p is an improvement over the trivial approximate 
solution. It is well known that there exist quasimodes for the linear equation localized near F of the form 

Vhix) = /i-(''-i)/4e'^/'7(s, h-^/^x), (0 < /i< 1), 
with / a function rapidly decaying away from F, and s a parametrization around F, so that satisfies 

{X + Ag)vh=0{h^)\\vh\\ 

in any seminorm, see jRal82j . Then 

(A + Ag)vh = (j\vh\Pvh + Q2{vh), 
where the error Q2{vx) = satisfies 

3. A TOY MODEL 

In this section we consider a toy model in two dimensions, and we give an idea of the proof of Theorem 
11.31 For simplicity, we moreover assume that 2 < p < 4. As it is a toy model, we will not dwell on error 
analysis, and instead make Taylor approximations at will without remarking on the error terms. Consider 
the manifold 

M = R^/2t:Z X Me /27rZ, 

equipped with a metric of the form 

where A € C°° is a smooth function, A> e > for some e. From this metric, we get the volume form 

dVol = A{x)dxde, 

and the Laplace-Beltrami operator acting on 0-forms 

Agf ^ {dl + A-^d^g + A-'A'd,)f. 

We observe that we can conjugate Ag by an isometry of metric spaces and separate variables so that 
spectral analysis of Ag is equivalent to a one-variable semiclassical problem with potential. That is, let 
S : L'^{X, dVol) ^ L'^iX, dxdO) be the isometry given by 

Su{x,e) = A^/^{x)u{x,9). 

Then A = SAS~^ is essentially self-adjoint on L^{X,dxd9) with mild assumptions on A. A simple 
calculation gives 

-Af = i-dl - A-^{x)d'g + V^{x))f, 

where the potential 

V,{x) = \A"A-'^l{ArA-\ 

We are interested in the nonlinear Schrodinger equation (|l.ip . so we make a separated Ansatz: 

ux{t,x,9) = e-'*^e*'^V(a;), 

where fc £ Z and is to be determined (depending on both A and k). Applying the Schrodinger operator 
(with A replacing A) to u\ yields the equation 

{Dt + A)e'*^e*'=V(2;) = {X + dl- k^A-^{x) + Vi{x))e''^e'''^ij{x) = a|V'|Pe'*^e''=V(^), 
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where we have used the standard notation D = —id. We are interested in the behaviour of a solution or 
approximate solution near an elliptic periodic geodesic, which occurs at a maximum of the function A. For 
simplicity, let 

A{x) = v/(l + cos2(a;))/2, 

so that in a neighbourhood of x = 0, ~ 1 — and A^"^ ~ 1 + . The function Vi{x) ~ const, in 
a neighbourhood of a; = 0, so we will neglect Vi. If we assume "0(2^) is localized near a; = 0, we get the 
stationary reduced equation 

Let h = \k\~'^ and use the rescaling operator Tip{x) — Thfiipi^x) — h~^^^^{h~^/'^x) (see Lemma [4.51 below 
with n = 1) to conjugate: 

T-\-X + dl- k^il + x^))TT-^^j = r-i((T|Vr^) 

or 

(-A + h~^dl - k\l + hx^)ip = ah~P^^\ip\Pip, 
where ip = T^^Tp. Let us now multiply by h: 

i-dl + x^-E)^^ah'>\ip\Pip, 

where 

h 

and 

1 d-1 p 

(7 := 1 — p = i . 

^ ^4 4 

Observe the range restriction on p is precisely so that 

0<q< 1/2. 

We make a WKB type Ansatz, although in practice we will only take two terms (more is possible if the 
nonlinearity is smooth as in the context of Theorems 11.11 and II. 2p : 

if = tpQ + h'^ifi, E ^ Eo + h'^Ei. 

The first two equations are 

: {-dl + x^ - Eo)'Po ^ 0, 
/i« : i~dl + x^--Eo- h'^Ei)^i = Ei^o + (t|(^o|Vo. 
Observe we have included the h'^Ei(pi term on the left hand side. 
The first equation is easy: 

^o(x) = e-^'/^ Eo = l. 

For the second equation, we want to project away from ipQ which is in the kernel of the operator on the 
left hand side. That is, choose Ei satisfying 

{Eiipo + (j\ipo\Pifio,ipo) = 0, 

so that the right hand side is in ipQ C L^. Then since the spectrum of the one-dimensional harmonic 
oscillator is simple (and of the form {2m + 1), m £ Z), the operator {—d^ + x'^ — Eq — h'^Ei) is invertible 
on ipQ C with inverse bounded by (2 — Hence for /i > sufficiently small, wc can find (pi € 

satisfying the second equation above (here we have used that (po is Schwartz with bounded norms). 
Further, since ipo is Schwartz and strictly positive, so is |(/3or<y5o, so by propagation of singularities, ipi is 
also Schwartz. In particular, both ipQ and ipi are rapidly decaying away from x = 0. 
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Let now ip{x) = T{ipo{x) + h'^(pi{x)), and observe that by the above considerations, ip{x) is 0{h°°) in 
any seminorm, outside an h^/"^^^ neighbourhood of a: = 0. Let u = e^'^^e^'^^ip{x) so that ||w||L2(d2;de) ^ li 
and u is 0{h°°) outside an h^^^~^ neighbourhood of x = 0. Furthermore, u satisfies the equation (again 
neglecting smaller terms) 

{Dt + A)u = h-^ThT-\-\ + dl- k^A~^{x))TT~'^u 
= h-^T{ah'i\T-^u\PT~^u + 0{h\T-^u\P+^) 
= a\uY'u + Q, 
where Q satisfies the pointwise bound 

Q = 0{h''\u\P+^). 

We now let u be the actual solution to with the same initial profile: 

(A + A)u = ct|u|Pm 

Set Th = hP, then with the Strichartz estimates of Burq-Gcrard-Tzvctkov [BGT04j we prove (see Proposi- 
tion [821) that there exists > so that 

l|w - w||l~([o,t,J;L2(m)) < Ch", 

and therefore we can compute: 

\\u\\L'-ilO,T^];L^M\U^,/2-6)) < \\u\\l'^{IQ,T]-L^M\U^i/2-s)) + 11"" " "!U~([0,T];L2(m\(7,^i/2_,)) 

which gives the result. 

Remark 3.1. It is very important to point out that the sources of additional error in this heuristic expo- 
sition have been ignored, and indeed, to apply a similar idea in the general case, a microlocal reduction to 
a tubular neighbourhood of T in cotangent space is employed. The function ipo is no longer so simple, and 
the nonlinearity |</3o|^Vo 'is no longer necessarily smooth. Because of this, the semiclassical wavefront sets 
are no longer necessarily compact, so a cutoff in frequency results in a fixed loss. 

Remark 3.2. We remark that the Strichartz estimates from |BGT04| are sharp on the sphere for a 
particular Strichartz pair, but this is not necessarily true on a generic Riemannian manifold. See |BSS08j 
for a thorough discussion of this fact. 

4. Preliminaries 

4.1. Symbol calculus on manifolds. This section contains some basic definitions and results from semi- 
classical and microlocal analysis which we will be using throughout the paper. This is essentially standard, 
but we include it for completeness. The techniques presented have been established in multiple references, 
including but not limited to the previous works of the second author |Chr07|[Chrllj . Evans- Zworski |EZ07j . 
Guillemin-Sternberg [GSTTI IGSTO] . Hormander |Hor03l IHotOS] . Sjostrand- Zworski [SZ02], Taylor |Tay81| , 
and many more. 

To begin we present results from jEZOTj . Chapter 8 and Appendix E. Let X be a smooth, compact 
manifold. We will be operating on half-densities, 

u{x)\dx\^ e (x,n\^ , 
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with the informal change of variables formula 

u{x)\dx\^ = v{y)\dy\^, for y = k{x) ^ v{k{x))\k' {x)\^ = u{x), 

where |K'(a;)| has the canonical interpretation as the Jacobian | det(9At)|. By symbols on X we mean the 
set 

= {aeC°°(T*Xx (0,l],nj.;,): |a^afa(a;,e;/^)| <C7„0/i-™(O'=-|^l}. 

Essentially this is interpreted as saying that on each coordinate chart, Ua, oi X , a = Qa where is Ua G S'''™' 
in a standard symbol on R*^. There is a corresponding class of pseudodifferential operators f2^) 
acting on half-densities defined by the local formula (Weyl calculus) in ]R"3 

Ovi;:{a)uix) = J J a e^^^-y-i>/My)dydC 

We wiU occasionally use the shorthand notations a"" := Op'^{a) and A := Op^{a) when there is no 
ambiguity in doing so. 

We have the principal symbol map 

which gives the left inverse of Op ™ in the sense that 

(T/i o Op ,j : 6 ' ^6 /6 

is the natural projection. Acting on half-densities in the Weyl calculus, the principal symbol is actually 
well-defined in 5fc,™/5fc-2,m-2^ ^^lat is, up to 0{h'^) in h (see, for example |EZ07| . Appendix E). 

We will use the notion of semiclassical wave front sets for pseudodifferential operators on manifolds, see 

Hormander |H6r05j . |GS77j . If a € iS'^'™(T*A', f2|,.j^), we define the singular support or essential support 
for a: 

ess-supp ha C T*X \_\ S*X, 

where §*X ~ {T*X \ {0})/R+ is the cosphere bundle (quotient taken with respect to the usual multipli- 
cation in the fibers), and the union is disjoint. The ess-supp /jO is defined using complements: 

ess-supp ha 

C[{x,OeT*X:3e>0, (5^af a)(x', O = d(x, x') + 1^ - ^'1 < e} 

|JC{(x,e) e r*X\0 : 3e > 0, (9,"af a)(x', ^ = O{h^{0~n, 

d(x,x') + i/ie'i + ic/iei-e7ie'ii<e}/M+. 

We then define the wave front set of a pseudodifferential operator A e \I',j'™(Ar, flj^): 

WFh{A) ess-supp h(a), for A = Ophi^)- 



^We use the semiclassical, or rescaled unitary Fourier transform throughout: 
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If u{h) is a family of distributional half-densities, u G C°°((0, 1]^, '^'{X, ^x)), we say u{h) is h-tempered 

if there is an A''o so that h^"u is bounded in ^'{X, ^x). If m = u{h) is an /i-tcmpcrcd family of distributions, 
we can define the semiclassical wave front set of u, again by complement: 

C{(x,0 : 3A e ^'If, with MA){x,0 + 0, 

and Au e h°°C°°iiO,l]hX°°{X,nl))}. 

For A = Op)^(a) and B = Op]^{b), a e gk.m^ b E S'^ wc have the composition formula (see, for 
example, [dMq] ) 

(4.1) v4oi? = Op^'(a#6), 
where 

(4.2) ^fc+fe'.m+m' 3 a^5(a;^^) ^^u{Dx.d,-d„d,) (^a{x , ^)biy , 7])) , 

with u} the standard symplectic form. 
We record some useful Lemmas. 

Lemma 4.1. Suppose (a;o,Co) ^ WF/i(w). T/ien V6 e C^(T*M") uii</i support sufficiently close to (a;o,^o) 
we /laue 

^(a;,/!!?)^ = 05(/i°°)- 

Here Os{h°°) means 0{h°°) in any Schwartz semi-norm. The proof of this Lemma follows similarly to 
that of Theorem 8.9 in jEZ07| . 

Theorem 4.2. (i) Suppose a G S{m) and u{h) is h-tempered. Then 

WFk{a'"u) C WFft(it) n ess-supph{a). 
(ii) If a Cz S{m) is real-valued, then also 



WFhiu) C WF,,(a"'ii) U C{ess-suppha}. 

Proof. Assertion 1 is straightforward. The proof of assertion 2 is standard, however we present it here so 
we can use it for the analogous result for the blown- up wavefront set. 

We will show if a{xo,^o) 7^ and a'^u = 0^,2 (/i°°) then there exists 6, 6(a;o,^o) 7^ so that b'^u = 
0]^2{h°°). There exists a neighbourhood U 3 (a;o,^o)i a real-valued function Xj and a positive number 
7 > such that supp x H C/ = and 

|fi + *xl ^ 7 everywhere. 
Then P = a"' -\- ix^ has an approximate left inverse c™ so that 

c'"P = id + i?"', 

where i?""' = OL^^L^-ih^). Choose b e S so that supp (5) C U and b{xo,^o) + 0. Then h^x"" = C(/i°°) as 
an operator on I? . Hence 

b'^u =b'"c'"Pu - b'"R'"u 

=b'^c'^a"'u + ib"'c'"x"'u - b'"R"'u 
=0(/i°°), 

where we have used the Weyl composition formula to conclude ess-supp ft (c#x) H supp 6 = 0. □ 
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4.2. Exotic symbol calculi. Following ideas from |SZ02j . since rescaling often means dealing with sym- 
bols with bad decay properties, we introduce weighted wave front sets as well. Let us first recall the 
non-classical symbol classes: 

= {aeC°°(T*Xx (0,1], r!|;.;^) : \d^d^a{x,e,h)\ < C„/3/i-*I"I-tI^I-™(0''"'^' } . 

That is, symbols which lose 6 powers of h upon differentiation in x and 7 powers of h upon differentiation in 
^. Note the simplest way to achieve this is to take a symbol a{x, ^) £ S and rescale {x, ^) 1— > {h^^x, h^^£,), 
which then localizes on a scale h^'^'^ in phase space. We thus make the restriction that < S,j < 1, 
< i5 -|- 7 < 1, and to gain powers of h by integrations by parts, we usually also require (5 + 7 < 1. We can 
define wavefront sets using Ss^-y = Sg'^ symbols, but the localization of the wavefront sets is stronger. 

Definition 4.3. // u{h) is h-tempered and 0<(5, 7<1,0<(5 + 7<1, 

for some a G 5 and d^u = 0^2 (/i°°)}. 

We have the following immediate corollary. 
Corollary 4.4. If a £ Ss^-y{m), < (5 + 7 < 1, and a is real-valued, then 



WFh^s,-yiu) C WFfi^s.-y{a"'u) U C{ess-supp ha}. 



The proof is exactly the same as in Theorem 14.21 only all symbols must scale the same, so they must 
be in Ss^-y- In order to conclude the existence of approximate inverses, we need the restriction S -\- j < 1, 
and the rescaling operators from ^ 14. 31 which can be used to reduce to the familiar h^^ calculus, where 
S' = {S + 7)/2, by replacing h with h^^'^ . 

4.3. Rescaling operators. We would like to introduce /i-dependent rescaling operators. The rescaling 
operators should be unitary with respect to natural Schrodinger energy norms, namely the homogeneous 
H" spaces. Let us recall in R", the space is defined as the completion of S with respect to the topology 
induced by the inner product 



where as usual u denotes the Fourier transform. The H'^ norm is just the norm, and the norm is 

The purpose in taking the homogeneous norms instead of the usual Sobolev norms is to make the rescaling 
operators in the next Lemma unitary. 

Lemma 4.5. For any s G R, /i > 0, the linear operator T^.s defined by 

Th,sw{x) = h'/^-^'/^wik'^/^x) 
is unitary on H^{M"), and for any r G M, 

Moreover, for any pseudodifferential operator P{x, D) in the Weyl calculus, 

Tf^lP{x,D)Th,s = P{h^/^x,h~^/^D). 
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Remark 4.1. Observe that for this lemma, the usual assumption that h be small is not necessary. 

Proof. The proof is simple rescaling, but we include it here for the convenience of the reader. To check 
unitarity, we just change variables: 

To check the conjugation property, we again compute 

= T-li2nr^h^/'+-/^ J e^{^-^'''y'^)p C^^^"\ 0^{y)dydi 

= T-l{2n)-h^f^-/^ J e<^~'^'''y''^^'''^)p{^±^,h-'^'OMy)dyd^ 

= {2n)-"l e^i^-y'i^P{l!—i^±yl,h-^f^^)n:{y)dyd^ 
^ P{h}'^x,h-^/^D)w{x). 

□ 

The purpose of using the rescaling operators Th^s is that if w S has /i-wavefront set 

WFft,5,7(") C {\x\ < a{h) and |^| < /3(/i)}, 

where, according to the uncertainty principle, a{h)/3{h) > h, then Th^sU has ft,-wavefront set 

^F,,^s-i/2n+i/2iTh^su) C {\x\ < a{h)h^'^ and |^| < mh'^'^). 

provided, of course, that 5 > 1/2. To see this, we just observe that for any %[) G C^(R^"), ip = on 
{\x\ < 1, ICI < 1}, we have 

^(x/a{h),D/m)u = Oih°°)\\u\\L2, 
or any other semi-norm, and hence 

= ThMx/a{h),D/l3{h))u 

^h''/^0{h°^)\\u\\L2 
^0{h°°)\\u\\L2. 

Finally, we note that the symbol of the operator on the left-hand side is zero on the set 

{|a;| < h^'^a{h), and |^| < h-^/^P{h)}, 



NONLINEAR QUASIMODES 



11 



and any such symbol in 55_i/2, 7+1/21 elliptic at a point outside this set, can be locally obtained in this 
fashion. 

5. Geometry near an elliptic geodesic 

Suppose r is a geodesic in a n + 1 dimensional Riemannian manifold. Following |MP05[ §2.1] (cf. 
jGra90| ) we fix an arclength parametrization ^{t) of F, and a parallel orthonormal frame Ei, . . . ,En for 
the normal bundle A^F to F in M. This determines a coordinate system 

X = (xo,xi, . . . ,a;„) 1-^ exp^(^^) (xiE'i + . . .x„£'„) = F{x). 
We write x' — (xi, . . . , x„) and use indices j, k,£ € {!,..., n}, a, /3, (5 G {0, . . . , n}. We also use = 



Note that r{x) = \Jx\^ . . .x\ is the geodesic distance from x to F and dr is the unit normal to the 
geodesic tubes {x : d{x,T) = est}. Let p = F{xo,0), q = F(xo,x'), and r = r{q) = d(ji,q) then we have 
[MP05[ Proposition 2.1], which states 

1 q 

gjkiq) = Sij + -g{R(Xs,Xj)Xi,Xk)pXsXi + ©(r^). 



(5.1) 



500(g) - l-g{R{Xk,Xo)Xo,Xe)pXkXi + 0{r^), 

71 

Too = -J29{R{Xk,X,)X,,X,)pX,+0{r^), 
i=i 



where F*^ = ig*'''(Xag^^ + Xpgar, - X^gap). 

In these coordinates, the Laplacian has the form 

1 

Vdetg 



Ag ^==div(v/d^5~^V) 



= .g°°^o^o + 2,9'=°XfeXo - .g^''=F°fcXo + 2.g°^FS,Xfc + Ap^ , 

where Apj_ is the Laplacian in the directions transverse to F. 

Denote the geodesic flow on SM , the unit tangent bundle, by (ft- Let 7(0) = p € F and C = 7'(0) S SM . 
Associated to F is a periodic orbit iptC of the geodesic flow on SM . This orbit LptC is called stable if, whenever 
V is a tubular neighbourhood of LptC,^ there is a neighbourhood U oi such that (' £ U implies (ptC ^ V. 
Given a hypersurface S in SM containing C and transverse to (^(C, we can deflne a Poincare map V near 
C, by assigning to each C,' the next point on (fitiC) that lies in E. Two Poincare maps of ftC are locally 
conjugate and hence the eigenvalues of dV at C are invariants of the periodic orbit LptC- 

The Levi-Civita connection allows us to identify TTM with the sum of a horizontal space and a vertical 
space, each of which can be identified with TM . Thus we can choose S so that TqYj is equal to E Q) E, 
where E is the orthogonal complement of C in TpM . The linearized Poincare map is then given by 

E®E — ^E®E, 

(y, W) I (J(length(7)), Vxo J(length(7))) 



12 



P. ALBIN, H. CHRISTIANSON, J.L. MARZUOLA, AND L. THOMANN 



where J is the unique Jacobi field along 7 with J(0) = V and Vx,, J(0) ~ W, i.e., J solves 

Nxo'^x„J + RiJ,Xo)Xo^O 
\j{0) = V, VxoJiO) = W. 

The linearized map dV — dVl^- preserves the symplectic form on E Q) E, 

Uji(yi,Wi), iV2,W2)) - {VuW2) - {Wi,V2), 

and so its eigenvalues come in complex conjugate pairs. We say that F is a non- degenerate elliptic closed 
geodesic if the eigenvalues of dV have the form {e^'"j : j = 1, . . . , n} where each aj is real and the set 
{oil, ■ • ■ , Q!„, tt} is linearly independent over Q. 

From (|5.1|) , the Hessian of the function goo (?) £is a function of x' is (minus) the transformation appearing 
in the Jacobi operator 

E3V ^ R{V, Xo)Xo e E. 
Notice that if y is a normalized eigenvector for this operator, with eigenvalue A, then 

seCp(Xo, V) = gp{R{V, Xo)Xo, V) = \gp{V, V) = A. 

The very useful property 

(5.3) peT,V eE =^ seCp(Xo, V) > 0, 

holds for any elliptic closed geodesic. 

Remark 5.1. One can verify (|5.3p hy means of the Birkhoff normal form (see [Zel08| §10.3]J. Indeed, 
if one of the sectional curvatures were negative, then the Birkhoff normal form of the linearized Poincare 
map (in T*M) must have an eigenvalue off the unit circle. If so, then there is at least one nearby orbit 
which does not stay nearby, hence the periodic geodesic is not elliptic. Similarly, if one of the curvatures 
vanishes, then the linearized Poincare map has a zero eigenvalue, and hence the logs of the eigenvalues are 
not independent from n over the rationals (in other words, the Poincare map is degenerate, and not even 
symplectic) . 



6. Compact Solitons: The nonlinear Ansatz 
We are interested in finding quasimodes for the non-linear Schrodinger equation 

— AgU = Xu — (t|u|^u, 

where a — ±1 determines if we are in the focussing or defocussing case. We will construct u approximately 
solving this equation with u concentrated near F in a sense to be made precise below. 

We take as Ansatz 

F{x,h) = e^^"/'^f{x,h), h-'^^, 

where L > is the period of F, and assume for the time being that the function / is concentrated in an 
/i-dependent neighbourhood of F. We are going to employ a semiclassical reduction, and we are interested 
in fast osciUations {h 0), so we assume WF^ 1/2-15,0/(2^; ^t-) C {\x'\ < eh^/'^~^ , |^'| < e} for some e, (5 > C0. 



■^Thc reason for the weaker concentration in frequency | is that the nonhnearity forces working with non-smooth functions, 
so some decay at infinity in frequency is lost. 
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The localization property of / as constructed later will be verified in Section [T] We compute from (|5.2p , 
with A the non-positive Laplacian, 



(6.1) AF 



co/'i 



g 



1 



2i 



+ jXof + XoXof ) + 2g'^°Xk ( -/ + Xof 



-.9 ^r' 



{h 



Remark 6.1. One may initially be inclined to use the Ansatz of the original Gaussian beam from Ralston 
[Ral82j ■ which is 

e'^('=y^{ao + aih + ---+aNh^), 

the standard geometric optics quasimode construction. After all, Ralston is able to make very nice use 
of the geometry to construct a phase function of the form i/h{xo + ^x' B{xo)x') with Im_B(xo) > (for 
Xq 7^ 0, vanishing otherwise). In such a regime, however, the non-linear term in the Schrddinger equation 
(|l.ip vanishes to infinite order in h. Thus while such a solution always exists, it fails to capture the effects 
of the nonlinearity that we are interested in. 

We analyze (|6.1|) by applying the operator Tf^^ in the variables transversal to F. We normalize every- 



thing in the L sense, so we take here s — 0. Let 



h ^^^x' and set v{xo,z,h) = Tj_^^f{xQ,z,h) 



h"'^^f{xo,h^^^z,h). Notice that the distance to the geodesic r = \x'\ is scaled to p = |z| = /i ^^^r, as 
described above. In particular, now 

(6.2) WF^,o,i/2^' C {|z| < h-'e, |C| < h'/h}, 

if C is the (semiclassical) Fourier dual to z. We conjugate (|6.ip to get 



^00 



2i 



/l2 



—V + -Xov + XoXov + 2g''%-^/^d,, ( -rv + X, 



0]' 



where the metric components and Christoffel symbols are evaluated at {xo,h^^^z). On the other hand, 
from (|5.ip . expanding in Taylor polynomials, we know that 

g™(a;) = 1 + Rsix) + Rsix) + Ri{x) + 0{r^) 

= 1 + hR2{z) + h^^^Rsiz) + h^Ri{z) + Oih'^^^p^) 

g'Hx)^h~gfiz)^0{h''^p') 

g=\x)^5,k + 0{hp') 

F°,(x) = h''^V%,{z) + hf^^^^iz) + 0(/.3/2p3) 
TUx) = h^/^Tl^,{z) + 0{hp^), 



for some smooth functions Rg, gi , homogeneous of degree £ respectively. Hence 



.3) T-lAn.oT-'.F 



Jxo/h 



1 



1 



1 



2i 



2i 



-l-g^'^r%,{z)v ~ ig^'T%,{z)v + ^Ar.v 



Pv, 



where P contains the remaining terms from the Taylor expansion. Let us record the following Lemma, 
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Lemma 6.1. The operator P has the following expansion: 

P = 0(/ll/2|z|5) + XoXo + 0(/l|zn/l-l/25^^Xo + 0(|zp)a,, 

In particular, if v satisfies (j6.2p , then 

\\Pv\\l- < C\\v\\l2 + C\\XoXov\\l2 + Ch'/^-'\\Xav\\L2. 

Remark 6.2. We will show later that for the particular choice of v we construct, the operators Xq and 
Xq are bounded operators, so that the error Pv is bounded in by v (see Remark \7.3\ ). 

For the purposes of exposition, let us then assume for now that the term Pv is bounded and proceed 
(this will be justified later). Applying Tj^^ to the equation — AF = XF — a\F\PF yields 

= XT-^F - ah'P-/^\T-lF\PT-lF, 
so that multiplying (|6.3p by /le"""/'', we get 

(6.4) {2iXo + Ar± ~ R2{z))v 

where P is the same as above. 

Remark 6.3. In order to ensure that the nonlinearity appears here as a lower order term, we require 

4 4 

(6.5) 9 := 1 —pn/i > 0, or p < — = — 

71 a — 1 

as stated in the theorems. 

We want to think of the left hand side as similar to a time-dependent harmonic oscillator where xo plays 
the role of the time variable. 

Let q = 1 — pn/A, < q < 1. We would like to assume that v and A have expansions in h'^, however 
the spreading of wavefront sets due to the nonlinearity allows us to only take the first two terms when 
< g < 1/2 and the first three terms otherwise. 

Case 1: < q < 1/2. Assume that v has a two-term expansion 

V = vq + h''vi 

and moreover that there exist Ek, fc = 0, 1, satisfying 

tlJf^=Eo + h'^E,+0{h''^). 

bmce q < 1/2, then the 0{h^/'^) term in (|6.4p is of equal or lesser order than the nonlinear term, and 
substituting into (j6.4p we get the following equations according to powers of h: 

h^ : {2iXo + Ar± R2{z) - Eo)vo ^ 0, 

(6.6) 
where 

Lvo = [tg'''f%, - 2i~g'^\z)d,, + R^{z) + h^/^R^{z) + ih^'^g"^^^^ + h^'^e~'-°"'p) vo- 
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We will show the error terms are 0{h?'^) in the appropriate H'' space. See ^7.2.11 

Case 2: 1/2 < g < 1. In the case q > 1/2, the 0{h^^^) term becomes potentially larger than the 
nonlinearity, so we take three terms in the expansions of v and E = ■ 

v = vo + /i^/^i'i + h'^V2, E = Eo + h^l'^Ei + h'iE2 + 0{h). 

We then want to solve 

(22Xo + Ar^ -i?2(-z)-£^o)wo = 0, 
h^l^ : {2iXo + Ar± - R2{z) - Eq)vi = EiVq + Lvq 
/i« : {2iXQ + A-p±-R2{z)-EQ-h^^^Ei-h''E2)v2^E2Vo + h^'iEiVi+h^^'^E2Vi 

+ a\vofvo + h^-''Lvi. 

In this case we will show the error is 0{h^^^^'^) in the appropriate H'^ space. See ij7.2.2l 

7. QUASIMODES 

We begin by approximately solving the h'^ equation by undoing our previous rescaling. That is, let 
wo{xo, X, h) = Th.o'Vo{xo,x,h) = h^'^/'^vo{xo, h~^^'^x, h), and conjugate the equation by Thfi to get: 

= ThfiC^iXo + Ar± - i?2(z) - Eo)T^}^ThfiVo 

^ {2iXq + /lAr^ - h-^R2ix) - Eo)wo, 

where now the coefficients in Api are independent of h, and we have used the homogeneity of i?2 in the x 
variables. Multiplying by h, we have the following equation: 

(7.1) (2ihXo + h^Ar^ - Raiz) - hEo)wo = 0, 

Hence, ()7.ip is a semiclassical equation in a fixed neighbourhood of T with symbols in the ft," calculus (i.e. 
no loss upon taking derivatives). The principal symbol of the operator in (|7.ip is 

p = r-\C\l-R2{z), 

where g is the metric in the transversal directions to F, and r is the dual variable to the xq direction. If 
we let r be the (unit speed) lift of F to T*M, and if exp(siJp) is the Hamiltonian flow of p, then 

r {c = z 0, .To = s e M/z}. 

Since, in the transversal directions, p is a negative definite quadratic form, the linearization of the Poincare 
map S is easy to compute: 

S = exp{Hq), 

where 

9=-|Cl^-^2(^), 

so that 

H, - -2a(To,z)^'^05., +26^''=(xo)z,%, 

where a and b are symmetric, positive definite matrices. Linearizing S about, say, xq and z = ^ = we get 
that dS{0, 0) has all eigenvalues on the unit circle, in complex conjugate pairs. That is, F is still a periodic 
elliptic orbit of the classical flow of p. 

Since p is defined on a fixed scale, we can glue p together with an operator which is elliptic at infinity 
so that p is of real principal type so that we can apply Theorem lA.ll in the appendix to construct linear 
quasimodes. Note that since we have quasi-eigenvalue of order 0{h), Theorem I A . 1 1 implies the quasimodcs 
are concentrated on a scale \z\ < h^^'^, |C| < h^^^. 
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This is made precise in the fohowing proposition. 
Proposition 7.1. There exists G L'^ , II^oIIl^ = 1, and Ea ~ 0(1) such that 

i2ihXQ + /i^Ar-L - R2{z) - hEo)wo = eo{h). 

Here the error eo(/i) = 0{h°°) E (or in any other seminorm), and wq has h-wavefront set sharply 
localized on T in the sense that if (p G ^1/2-5 1/2-5 ^ riear T, then for any S > 0, (pwo = wq + 0{h°°), 
and if 6 = 0, \\(pwo\\ > Co||woll for some positive Cq depending on the support of (p. 
Moreover, Wq G H°°{AI) and satisfies the estimate 

\\wo\\H~^iM)^Oih''/^). 

Proof. The construction foUows from Theorem lA.ll and is well-known in other settings, see for instance 
jRal82j ■ |Bab68j . and |CP02| . To get the sharp localization, apply Lemmas IB . 1 1 and IB .21 from the appendix 
to get the localization on wq. Once we know that Wq is so localized, we can replace Wq with ipwo, where 
ip € proposition. Then (pwo satisfies 

{2ihXo + /i^Api - R2{z) - hEn)'pwo = eo{h), 

where 

eoih) = (peo{h) + [{2ihXo + h^Ar± - i?2(z)), ^]«^o- 
But now ipeo{h) = 0{h°°), while the commutator is supported outside of an h^/'^~^ neighbourhood of F, 
so by the localization of is 0{h°°) and localized in a slightly larger set on the scale h^^^~^ . □ 

Now recalling vq = Tf^^wo, then vq satisfies 

{2tXo+Ar± - R2{z) - Eo)vo 

= h~^T-^{2ihXo + h^Ar± - R2{z) - Eo)nflVo 
= h-^T-l{2ihXa + h^Ar± - i?2(z) - Eq)wo 
^ h-'T^',eo{h). 

The error /i^^r^^eo(/i) = 0{h°°) in any seminorm still, but the function vq is now localized on a scale 
in space. That is, 

WFh^o,i-svo C {\x\ < eh-\ < eh^-^}. 

7.1. The inhomogeneous equation. We are now in a position to solve the lower order inhomogeneous 
equations in (|6.6p . The quasimode Vq has been constructed as a "Gaussian beam" (see |Ral82| ): it is a 
harmonic oscillator eigenfunction extended in the direction by the quantum monodromy operator from 
[SZ02| . which is defined in (|7.4p below. From this construction, the boundedness of the error term Pvq 
as stated in Lemma 16.11 In what follows we construct vi in the case < g < ^ and vi , V2 in the case 
i < g < 1 and show that a similar bound holds for h'^Pvi and h^Pvi + h'^Pv2 respectively. 
We want to solve 

(7.2) {2iXo + Ar-L - i?2(z) - Eq)vi = Eivq + G, 

where Eq ~ Eq + h'^Ei for some 7y > 0, i?o and vq have been fixed by the previous construction, and 
G ~ G{xo,x) is a given function (periodic in xq) with 

WFh,o,i/2Gc{\x\<eh-',\^\<eh'/^}. 

Note, the localization of G follows from the nonlinearity as well as geometric multipliers in the operator 
L, see Appendix |E] and (|6.2p . 
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Conjugating by Th,o as before, we get the equation 
(7.3) {2ihXo + /i^Ar-L - R2{x) - hEo)wi = Gs, 

where 

G2 = hTh,a{Eivo + G). 

We observe then that 

WF,,i/2_5,oG2 C {\x\ < eh^'^'\ lei < e}. 

Spccificany, let M{xo) be the deformation family to the quantum monodromy operator defined as the 
solution to: 

{2ihXn + /i2Ar_L - i?2(x) - hEa)M{xQ) = 0, 
M{xo) = id, 

which exists microlocally in a neighbourhood of F C T*X (for further discussion see also Appendix lAl and 
references therein). By the Duhamel formula, we write 

wi = M{xo)wi^o + T Af(xo)Af (yo)*G2(yo, O'^yo- 



(7.4) 



We have to choose wi.o and Ei (implicit in G2) in such a fashion to make wi periodic in xq; in other 
words, to solve the equation (approximately) around F. Let L be the primitive period of F, so that xq — 
corresponds to — L. Then we require 

wi{L,-) =wi{0,-), 

or 

wi,o=^ M{L)wi,o + ^ / Af(i)Af(yo)*G2(yo,-)dyo. 



h Jo 

In other words, we want to be able to invert the operator (1 — M{L)). The problem is that wq.o •= 
u>o(0, •) = ThfiVo{0, •) is in the kernel of (1 — M{L)), so we need to choose Ei in such a fashion to kill the 
contribution of G2 in the direction of wo,o- 

Recall that 

G2 ^ hThfliEiVo + G) 
= h{Eiwo + G), 

where 

G = Tfi.oG. 

We want to solve microlocally 

(1 - M{L))wifl = - / M(L)M{yQy{h{EiWo + G))dyo 

L 

M{L)M{yonE,M{yo)wo.o + G)dyo 
(7.5) ^iLM{L)EiWofi+i f M{L)M{yoy Gdy^. 



Let 



1 / 



Ei^--{ I M{yoyGdyo,wo,o 



so that (|7.5p is orthogonal to wo,o • 
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If we denote 



= {u e : (it, wo,o> = 0}, 



then by the nonresonancc assumption (since Eq + li'^Ei is a smah perturbation of Eq), and the fact that 
M{L) is unitary on , {I — M{L))^^ is a bounded operator (see jChrllj ) 



{l-M{L))-^:Ll 



Hence 



satisfies 



M(L)M(yo)*G2(2;o,-)rfyo 



\\wi.o\\L^<Ch-' 



and 



Furthermore, we have the estimate 



M{yoyG2{yo,-)dyo 



Win e L 



2\\L^(xa)L'^ 



<CL^'^h-^-''HG 



L2 



that is, the norm is controlled, but not by the homogeneous Sobolev norm. 

We have proved the following Proposition, which follows simply from tracing back the definitions. 

Proposition 7.2. Let vq be as constructed in the previous section, and let G £ for s > sufficiently 
large satisfy 

WF„,o.i/2G C {|x| < e/i-MCI < 
Then for any j] > 0, there exists vi G and Ei ~ 0{\\G\\]^2) such that 

i2iXo + Ar-L - R2{z) - Eq - h'^Ei)vi ^ E^vq + G, 



ana moreover 



\\vi\\h^<C{\\G\\h^+\\vo\\h^). 

Remark 7.3. We note here that by construction of vi we have implicitly microlocalized into a periodic 
tube in the xq variable. Using the fact that the Quantum Monodromy Operator is a microlocally unitary 
operator (see jSZ021 [Chrllj ]. the bound 

\\X^viU.<\\vi\\r^. 
follows easily for any j , which is required for proof of Lemma \6.1i 

7.2. Construction of quasimodes in the context of Theorem 11.31 We now have all the tools to 
construct the quasimodes which will be used to prove Theorem 11.31 Let d > 2 and 0<p<4/(d— 1). As 
previously, denote hy q = 1 — p{d — l)/4. The main results of this part will be stated in Propositions 17. l1 
and [721 
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7.2.1. The case < q < 1/2, i.e. < p < In this subsection, we see how to apply Proposition 17.21 

in the case < q < 1/2. As described previously, in this case, the nonlinearity is the next largest term, 
and we have only one inhomogeneous equation so solve (see (|6.6p ). 

According to Propostion IE . 2 1 and Corollary IE. 31 from Appendix |E1 if 

is the nonlinear term on the right-hand side, then G is sharply localized in space but weakly localized in 
frequency. That is, if x C^{T*X) is equal to 1 in a neighbourhood of F, then for any < <5 < 1/2 and 
any < 7 < 1, 

Xih^x, h'">D,)Gixo,x) ^ Gixo,x) + E, 

where for any < r < 3/2, 

We are interested in the case where 7 = 1/2, since in that case G is weakly concentrated in frequencies 
comparable to h^^"^, so by cutting off, satisfies the assumptions of Proposition 17.21 That is, take 7 = 1/2, 
and replace G with G = x{h^ x, h^^^ Dx)G, and apply Proposition 1 7 . 2 1 to get vi and Ei satisfying 

{2iXo + Ar± - i?2(z) -Eo- h"Ei)vi = Ervo + G, 

or in other words 

{2iXo + Ar± - R2{z) - Eq - h''Ei)vi = E^vo + <j\vo\Pvo + h^^^-^Lvo + Qi, 

where 

IIQlllff. <C/l(l/2)(3/2+P--). 

Now letting v — vq + h''vi and E = Eq + h'^Ei, we have solved 

(2^X0 + Apx - R2iz) ~ E)v = /iVIwoTwo + h^l'^LvQ + /I'Qi 

where 

Q2 = h'^Qi - h^^^+^Lvi + 0(/i29|i,|P+i). 

The remainder Q2 satisfies 

since q < 1/2. Recalling the definitions, ip = e'^°^'^Th^ov satisfies 
(i) 

WF;,,i/2-M'P ^ {\^\ < ^h^'^'\ 1^1 < e}; 



(ii) 

and 

(hi) 



or 



= h-^e''''°/''Thfi{2iXo + Ar± - R2 - E - h^^^L)v 
= h-^e'^"^''n^oi<jh''\v\Pv + Q2), 

{A + X)^ = a\ip\Pip + h''^P^Q, 
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where Q = h-^e''-^°/''ThfiQ2 satisfies < Ch^' and where 

(7.6) a{p):=~l + 2q=l-p ^ ^ 



Wc now sum up what we have proven in a proposition. Consider the objects we have just defined : 

v^vo + h^vi, (ph = e"»/''r,i,ow and A(/i) = h-'^ - Eoh'^ - Eik'^+'i. Then we can state 

Proposition 7.1. Let 2/(d— 1) <J'<4/((i— 1) and a{p) he given by (j7.6l) . Then the function (ph satisfies 
the equation 

(A + \{h))^n = (T\vh\''Vh + h'^^P^Qih), 
where Q{h) is an error term which satisfies ^ ^ '^' for all s > 0. 

7.2.2. The case 1/2 < q < 1, i.e. < p < We again construct Vq as a Gaussian beam using the 

quantum monodromy operator. We then set 

Gi = Eivo + Lvo, 

which is smooth with compact wavefront set contained in the wavefront set of vq, so no phase space cutoff 
is necessary to apply the inhomogeneous argument to get vi with wavefront set contained in the wavefront 

set of Vq. 

Now let 

Ga = £^2^0 + h^-'^EiVi + h^'^E2Vi + a\vo\Pvn + h^-'^Lvi, 

and solve for V2 as in the previous subsection. This is possible since vi is orthogonal to vo by construction. 
We then have 

{2iXo + Ar_L - R2{z) - Ea - h'^Ei)v2 = E2V0 + h^-'^E^vi + h^/^E2Vi + a\vo\Pvo + h^~'^Lvi + Qi, 
where 

Letting v ^ vq + h^/'^vi + h''v2 and E = E^ + h^/'^Ei + h''E2, we have solved 

{2iXo + Ar^ - R2{z) - E)v ^ /iVIwoTwo + h^''^LvQ + HLvi + Qi 

= h'^a\v\Pv + h^/^Lv^Q2, 

where 

We now have the remainder estimate 

Recalling the definitions, ip := e^^°^^ThfiV satisfies 
(i) 

WF^^i/s.^-^oV' C {|x| < < e}; 

(ii) 

and 
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(iii) 

(A + A)^ = /i-ie"°/''T„,>e-"°/''T^-;>T^„r-> 

or 

(A + AV = a|</j|P(^ + /i"(p)Q, 
where Q = /i~^e""/''T/,,Q2 satisfies i|Qj|^, < Ch'" and where 

(7.7) a{p):=-ll2 + q=\-p{^y 

Once again, by construction we have 

for j = 0,1,2. 

Consider u uq + h^vi + /i«W2, i/S/i = e""/''T,j,of and A(/i) = h-'^ - -Bq/i^^ - - E2h-'^+i defined 

previously, then we have proven 

Proposition 7.2. Let Q < p < 2/(d— 1) a^d a{p) he given by (j7.7p . T/ien t/ie function iph satisfies the 
equation 

(A + A(/i))^;, = + /^"''''Ql/i), 

where Q{h) is an error term which satisfies \\Q{h)\\jj,, < h^", for all s > 0. 

7.3. Higher order expansion and proof of Thereorem II. ll 

7.3.1. The two dimensional cubic equation. We first deal with the simpler case d = 2, p = 2 and s = 0. 
As in the previous section, we define q = 1 — p{d — l)/4, thus for this choice we have q = \ allowing us 
to match powers of the asymptotic parameters in a canonical way. The general algorithm for any smooth 
nonlinearity arising when rescaling in the appropriate H'^ space will follow similarly. 

Using (|6.4p and the Taylor expansions of the geometric components g"^^ and r*^, for z, j', fc = 0, . . . , d in 
(|6.3p . we look for an asymptotic series solution of the form 

U = UO + h^Vi + h^V2 H h h'^^rn H 1- h^^VN + V, 

for N sufficiently large. 



22 



P. ALBIN, H. CHRISTIANSON, J.L. MARZUOLA, AND L. THOMANN 



Then, we have the following equations: 



{2tX„ + - R2iz) - Eo)vo = 0, 

{2iXo + Ar± - R2{z) - Eo)vi = Eiv^ + {i5jkf%i ~ ^w''"d,, + R3{z))vo + <j\vo\^vo, 
{2iXo + Ar± - R2{z) - Eo)v2 = ^2^0 + E^vi + {i6jkf%, - 2ig''°d,, + R3{z))vi 



/i" : {2iXo + Ar± - R2{z) - Eo)v^ = ^ E,n~jVj + cr( X! i^Tki^J^kVi)) 

m—1 
j=0 



where 



N N N-1 

{2iXo + Ar± - R2{z) - Eq^Y^ hi Ej)vn = E ^^-J^J + E i^jkiVjVkVi)) 

j=l j=0 j,k,l=0 

N-l 



J j.m 



ONi\zr-'), 
ONi\zr-^), 
fl,,, ^ ONi\zr-n 

for J, TO = 0, . . . , iV and 

= 0{h''/'-'\z\^) + XoXo + Oih^/'\z\^)h~'/'d.,Xo 

+0(;jW/2-3/2|^|iV)5^^ + + 0(/i^/2|z|^)Xo + 0(/i^/2|z|^)/i-i/2a,^_. 

Note that all constants have implicit dependence upon N relating the number of terms in the expansion 
at each order. The expansion is valid provided first of all that 

N 

Y,hiE,<Eo 

in order to justify the solvability of the 0{h^^'^) equation. 

Remark 7.4. We note here that in this expansion, the sign of a can effect the sign and value of Ei, 
which will impact the remaining asymptotic expansion and in particular the order of quasimode expansion 
possible. It is possible that the focussing /defocussing problem enters in to the stability analysis of these 
quasimodes through this point. 



Applying Proposition [72] at each asymptotic order and bounds similar to those in Lemma 16.11 at order 
as in Section 17.2.11 we have by a simple calculation that u is a quasimode for the nonlinear elliptic 
equation with remainder Qn such that 

\\QN\\H^<CNh~'-^+''l\ 
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As a result, for sufBcicntly smooth nonlinearities, one is capable of constructing higher order asymptotic 
expansions and hence a quasimode of higher order accuracy. 

7.3.2. The general case. Let d > 2, p € 2N and s > 0. Wc define here f/s = 1 + p{s — ^j^). Assume that 
p{^^ — s) < 1, or equivalently that > 0. Firstly, write w = h'^v. Then v has to satisfy (|6.4p but where 
the power in front of the nonlinearity is . Hence, we can look for v and E of the form 

N N 

v= h^'^'"^'^" Vj,e + v and E ^ h^/^+^i^ Ej^e + E. 

Since qs > 0, the nonlinearity does not affect the equation giving vq,q, and we have 

i2iXo + Ar± - R2{z) - £;o,o)wo,o - 0, 

which is the same equation as before. Then, using Taylor expansions, we write all the equations, similarly 
to the previous case, in powers of h^/^+^i^ . Again, we can solve each equation and obtain bounds of the 
solutions and of the error terms. Moreover, it is clear that we can go as far as wc want in the asymptotics, 
so that we can construct a 0{h°°) quasimode, and this proves Theorem ll.il 



8. Error estimates 



8.1. The regular case. In this section, we assume that p is an even integer. 



Fix an integer k > d/2 (the fact that k is an integer is not necessary). We then define the scmiclassical 
norm 

11/11^. = 11 (1-/^2 A) ''^/VIIl^(m). 

In the previous section we have shown the following : Given a e M, there exist two functions (ph S H^{M) 
and Q{h) € H^{M) and \[h) g R so that 

(A + X{h))^h + o\^>h?^h + h'^Qih). 
Moreover, microlocally the function takes the form 

(8.1) (^,,((7, x',h) = h-^+'e'^/''fia, /j-i/V, h), 
and we have \\Q{h)\\ij'' < C. 

h 

We set Uapp(t, •) = e~**'^^''Vft- Then if we denote by Q{h) := e^^*^Q{h), the following equation is satisfied 

(8.2) «9tMapp - AWapp = 0-|Uapp|^Mapp + h°'Q{h). 

Proposition 8.1. Let s > 0. Consider the function tph given by (|8.ip . Let u he the solution of 

I idtu — Au = (jImI^u, 
I u(0, •) = Ph- 



(8.3) 



Assume that a > — i h s + p{ h s). Then there exists C > and cq > independent of h so that 

Uapp||ioc([o,T,j;H=(A/)) < C/^(''+l)/^ 

forQ<Th<CQ hP(^-'') ln(i). 

This result shows that Uapp is a good approximation of u, provided that the quasimode ph has been 
computed at a sufficient order a. 
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Proof. Here we follow the main lines of |Tlio08a[ Corollary 3.3]. With the Leibniz rule and interpolation 
we check that for all / e H''{M) and g S W'''°°{M) 

(8.4) ll./,9llff,'; < ll/ilH,^!l,g!lL=^(M) + !l./llL^(M)ll(l-^'A)'/^g|U.(^^^^ 
Moreover, as fc > d/2, for all A, /a £ H''{M) 

(8.5) \\flf2\\H-^<h-''/'\\h\\H^Jh\\H^^. 

Let u be the solution of (|8.3p and define w ^ u — Uapp. Then, by (|8.2p . w satisfies 

idtW - Aw = a{\w + Uappr(w + Wapp) - lUappl^Mapp) - h°'Q{h) 

w{0,x) = 0. 

We expand the r.h.s. of (|8.6p . apply the operator (l — /i^A)*"^^ to the equation, and take the L^- scalar 

l. In 

product with (l - /i^A) ' w. Then we obtain 

P+i 

\in\\ ttl- ^ 
dt' 



(8.6) 



(8.7) < W^' <p''\\h}; + /i". 



We now have to estimate the terms \\w^ w^dd "'IIh''' ^o'' 1 — J !i P + L From (|8.4p we deduce 

(8.8) <||u;^||^.||<+i-^|U^(M) + |k^^ 
By (|8.5p . and as we have 

(8.9) II<p+p^"^IIl==(m) < 11(1 - h'Af\i;^-m,^^M) < 

thus inequality (|8.8p yields 
Therefore, from ()8.7p we have 

Observe that ||u)(0)j| tft = 0. Now, for times t so that 

h 

(8.10) h-^r>/2^\wr+i <,^p(-^+s)^^^^^ 

I.e. \\w we can remove the nonlinear term in (jS.Sp . and by the Gronwall Lemma, 

h 

(8.11) ||«;||^. < C/i^-'^t-^+^'e^''''"'^^^'*. 
If Co > is small enough, and t < cq/iP^^"^) In ^, then 



so that inequality (j8.10p is satisfied. By the usual bootstrap argument, we infer that for all 

Co/ 



h 



we have 

Finally, by interpolation we get Huill^f < '^ll^llif'' j hence the result. □ 
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We are now ready to complete the proof of Theorem II .21 Consider Uh, the exact sokition to (jl.ip with 
initial condition (ph, then by the previous proposition and the description of Uapp, we can write 

\\'^h\\L-°{[0,Tn];L^M\U^i/2-s)) ^ II ""app II ([0,T];I,2(M\C/^i/2-5 )) + W^h ~ Uapp\\ {[0,T];L^M\U^i/2-s)) 

(8.12) = 0(/i°°) + ©(/i^'^+i)/'*) = ©(/i^'^+i)/"), 

which was the claim. 

8.2. The non regular case and d = 2. 

In this section we compute the error estimate in the case of a non smooth nonlinearity in dimension 
d = 2. Moreover we restrict ourselves to the case s = in (|8.ip (case of an L^-normalized initial condition). 

Proposition 8.2. Let iph be the function given by (jS.ip with s = 0. Let u be solution of 

idtu — Am = (t|u|^m, 

u{0, •) = (fi. 

Let e > 0. For p G (0,4)\{1}, we set — h^, and in the case p = /i^^*^. Then there exists C > 

and V > independent of h so that 

||m - ■Uapp||Loo([o.Th]:L2(jV/f)) < Ck" . 

Remark 8.3. Note the difference between the results of Propositions (jS.ip (when s = and d — 2) and 
(|8.2p . In the first case, we have Th of order h^^^, which is better than h^ obtained in the second 

result. However, in this latter result, there is no restrictive condition on the size of the error term in the 
equation. 

Proof. First, we follow the strategy of Burq-Gerard-Tzvetkov |BGT041 Section 3.] Let < p < 4, choose 

r > max(p, 2) and take 1 < s < 1 (there will be an additional constraint on s in the sequel). Then 

r 

take q so that 7 + ^ = 5 and si = s — K For T > define the space 

Y" ^C{[0,T];H''{M)) n L''{[0,T];W'''{M)), 
which is endowed with the norm 

hlly. = max^|lM(t)|U» + ||(l-A)«i/2M|li.([o,j,j^i,). 

By the Sobolev embcddings, we have C L''([0,T],i°°). Now, define w — u~ Unrm. Thcu w satisfics the 
equation 

, . f idtW - Aw = a{\w + Us,pp\P{w + Uapp) - IWappTWapp) - h°'^P'^Q{h), 

\w{Q,x) = 0, 

with a{p) = 1 - p/2 when 2 < p < 4 and a{p) = 1/2 - p/A when < p < 2 (see ([73)1 and ([7?71) ) and 

||Q(/i)||h= < Ch-\ 

• Case < p < 1. In [CFHlOi Estimate (2.25)], Cazenave, Fang and Han prove that for all < s < 1 

(8.14) II |w + Uapp\^{w + Uapp) - \u app\^ U app\\ ^ , < C\\Uapp\\ H- + C\\w\\h- (||Wapp|lioo + ||w||^oc). 

Indeed, in |CFH10| . the estimate is not stated exactly with these indices, but the proof still holds true. 
Moreover, in |CFH10| . (|8.14p is proved for x G M'', but the inequality can be adapted to the case of a 
compact manifold thanks to a partition of unity argument. 
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Assume that w satisfies the equation 

idtw — Aw = F, w(0, x) ~ 0, 

then with the Strichartz estimates of jBGT04j . the estimate < C||F||^i/^s holds true. Thus, with 

p 

the notation 7 = 1 , with (|8.13|) and (|8.14p we have 

r 

My' < C r \\uapp\\H4M\l^+C r \\w\\H^{\\u,pp\\l^ + \\w\\l^)+CTh"^P^\\QU^H' 
Jo Jo 

< CT^\\Uapp\\L^H4M\Ys + C\\w\\Y^{T\\UapXL¥L^ + T"'\\w\\P..) + 

(8.15) < CT''h-'\\w\\'Ys + C\\w\\y^ {Th~P'^ + T4\w\\Py^) + CTh-'+^f'^~P/\ 
Similarly, we obtain 

(8.16) \\w\\l^l^- < CT^Wwry. + C\\wh^L^{T\\uapp\\l^L- + T^|lz«|l?-.) + CTh'/^~P^\ 
Therefore, if we define the semiclassical norm j| ||y^s by 

(8.17) \\u\\Y;^^h-^\\uh^L^ + \\u\\Y^, 
thanks to (|5J5)) and ((5?TB| we infer 

(8.18) \\w\\ys < CrT'/i-^||w||P^. +C||u;||Y,.(T/i-P/4 + rT'||«;||P^.) +(:7T/i-^+i/2-P/4. 

Next we use the inequality ab < — e^^a^^ H £-P2^p2 -^^i^ich holds for a,b,e > and 1 = 1. With 

Pi P2 Pi P2 

a suitable choice of e and pi (here we use that < p < 1) we get 

(8.19) CT-^h-'\\w\\l-s < ]:\\w\\Yi^ + C{T^h-'f'^^-P\ 
Now, re-inject (|8.19p into (|8.18p and obtain 

(8.20) \\w\\y^ < C\\w\\Y^{Th-P^^ + T'>\\w\\Py.) +CiT'fh~'y^^^-P^ +CTh-'+'^^-P/\ 

1 p 

We now perform a bootstrap argument : Fix e > and set ~ h^. Fix 1 < s < 1 . Then it is 

r r 

possible to pick > small enough so that 7 = 1 — - >sH — ^ — . Assume that 

r p 

(8.21) \\w\\y^ < h~^+\ 
Then 

which tends to with h, thanks to the assumption made on 7. Hence, for h > small enough, with (|8.20p 
we get 

Finally, observe that — s + 3p/4 + 1/2 > — ,s + v, and the assumption 7 > s H ^- — is equivalent 

P ^ 

to (p7 — s)/(l — p) > — s + V. Hence for /i > small enough, we recover HwHys < —h^'^^'^, and by 

^ 2 

the usual bootstrap argument, the condition (|8.2ip holds for Th = hP . Now we can deduce the bound 
< /i^ll^fllF^ < h'^ 1 which was the claim. 
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• Case 1 < p < 4. Here we have, by |CFH10[ Estimate (2.25)], for all < s < 1 

1 1 1'^ H" '^app 1^ (^W -\- Uapp) I'^^app l^^app 1 1 //s — 

CW,A\h^{\\u,,Xl^' + \\w\\l-2)\\wU^ + C\\wU^{\\u^,Xl^ + \\w\\l^). 
With the same arguments as for (|8.15p we get, with 7=1 — 1/?' 

\\w\\y- < CT^\\Uapp\\L^HA\Uapp\\^L^L--\\M\Y'^ + CT'^\\Uapp\\L^ hAM\y^ + 

+C(Th,pp||^^^^ +r^||w;||?..)||w;||y. +Cr/i"(P)||Q|U5,ff= 

with a{p) = 1 - p/2 when 2 < p < 4 and a{p) = 1/2 - p/4 when 1 < p < 2 (see (|7!6l) and (|7Jl) ). Then, by 
the same manner, we get the following a priori estimate with the semiclassical norm || (recall definition 

) 

(8.22) |ju;||y^= < C{T'<h-'-^'P-^'^''^ + ^-^'^1^^^ + CT'^h-'WwWl-s + CT^WwW^^ + CTh-'+°'^P\ 

We now perform the bootstrap : Let r > max (2,p) (there will be an additional constraint on r). Fix 1 < 

r 

s < 1 and set Th = hP. Then if r is large enough (recall that 7 = 1-1/r), the term T^^h'''~''P-'^^^^ +Thh-P^'^ 
tends to with h, therefore if ft, > is small enough, from (|8.22p we deduce that 

(8.23) \\w\\y^ < CT;ih-'\\w\\'y. + Cr^||w||^t' + CThh-'+'''^P\ 
Choose < < p + q;(p). As previously we assume that 

(8.24) \\w\\y,; < h-'+''. 
Then with IJjS^ we get 

II II ^ 

Next when r > is large enough (and under the assumption Q < v < p + a{p)), we have ^p~ s+p[—s + v) > 
~s + v, ^p+{p-\-l){—s + v) > —s + v and p— .s + a(p) > ~s + v. To see this, observe that 7 — )■ 1 and s — > 1 

when r — > +00. Therefore for /i > small enough, we recover < -h~'^^'^, hence the condition 

(|8.24p holds for Th = /i^, and similarly to the previous part, we deduce that ||M||iooi2 < /i" 1 1 w 1 1 < h"^ . 

• Case p = 1. By (|8.22p we have 

\\w\\y^ < C{T'-^h-' +Th-^/^)\\w\\Y^+CT'-r\\w\\l..+CTh-'+i. 
here we set Th = h^^'^ with e > 0, and we perform the same argument as in the previous case. □ 

Thanks to this proposition and the same argument as (j8.12p . we can conclude the proof of Theorem [LS] 

Appendix A. Quasimodes for linear equations near Elliptic Orbits 

In this section, we state, without proof, a theorem on existence of quasimodes near elliptic periodic 
orbits of the Hamiltonian flow. A proof can be found in jChrll] . 

Let X be a smooth, compact manifold, dimX — n, and suppose P £ 4'''-°(X), fc > 1, be a semiclassical 
pscudodifferential operator of real principal type which is semiclassically elliptic outside a compact subset 
of T*X. Let $i = cxptHp be the classical flow of p and assume there is a closed elliptic orbit 7 C {p = 0}. 
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That 7 is elliptic means if C {p = 0} is a Poincare section for 7 and S : N ^ S{N) is the Poincare map, 
then dS{0,0) has eigenvalues all of modulus 1. We will also need the following non-resonance assumption: 

if e=^*"S e=^'"^ . . . , 6^*"*= are eigenvalues of dS{Q, 0), then 
ai, a2, • ■ • , ttfc are independent over ttZ. 



(A.l) 



Under these assumptions, it is well known that there is a family of elliptic closed orbits C {p — z} for 
z near 0, with 70 — 7. In this work we consider the following eigenvalue problem for in a neighbourhood 
of z = 0: 



(A.2) 

We prove the following Theorem. 



{P-z)u = 0; 
\\u\\l^x) = 1- 



Theorem A.l. For each m G Z, 771 > 1, and each cq > sufficiently small, there is a finite, distinct 
family of values 

{z,}fj'l^C[-coh'^"\coh'^-] 
and a family of quasimodes {uj} ~ {uj{h)} with 

WFhUj = 7^^. , 

satisfying 

(A3) / {P - ^ 0{h°-)\\u,h^x); 

I \\uj\\mx) = 1- 

Further, for each m ^ "L, m > \, there is a constant C ~ C(co, 1/m) such that 

(A.4) < ^^(/i) < Ch-". 

Remark A.l. The proof is essentially to construct quasimodes on the Poincare section as eigenfunctions 
for the semiclassical harmonic oscillator, and then to propagate them around the orbit T. This shows that 
the quasimodes have the localization property as in Lemmas \B.l\ and \B.^ 



Appendix B. Some commutator estimates 



In this section, we prove two results which we have used in the above computations. Specifically, we 
have constructed approximate solutions to the homogeneous and inhomogeneous equations associated to a 
semiclassical operator of the form 

Q^hDt~P{t,x,hD,), 
where P is a "time-dependent" harmonic oscillator, 

Pit,x,hD,) = ~h^A{t) + V, 

where 

V = b'^{t)x,x, 

is a positive definite quadratic form. Both A{t) and V have i-periodic coefficients of period T (the same 
as F), and we seek periodic solutions to equations 

Qv = Ev, {Q~E)v^f, 

where E is an eigenvalue to be determined and / is periodic in t. The constant-coefficient semiclassical 
harmonic oscillator is well-known to have eigenfunctions of a semiclassically scaled Hermite polynomial 
times a semiclassical Gaussian. These eigenfunctions, for small eigenvalues, have semiclassical wavefront 
set at (0,0). Moreover, for any €,S > 0, if |(x,^)| > e/i^/^~*, these eigenfunctions are 0{h°^) in the 



NONLINEAR QUASIMODES 



29 



Schwartz space. The purpose of this section is to prove that for the periodic orbit case, similar localization 
occurs. 

Lemma B.l. Suppose v solves 

'Qv = Ev + 0{h°°)\\v\\, E^O{h), 
< 11^11 = 1, 
viO) = viT) ^ vo, 
where Vq satisfies the localization property: 

V(5,e>0,3VjeC,°°(M),?A = l on{\{x,^)\ < eh^/^-^}, 
with suppV' C < 2eh^/^^^}, we have 

OphWvo=vo + 0{h°°)\\vo\\. 

Then there exists Ci > 0, ei — > as e — > 0, such that for all < t < T, if ^ ^ C^(R), ^^{r) = 1 for 
{\r\ < 1} with supp'0(r) C < 2}, then 

Op^,{^{pit,x,0/{elh^-''))Ht)=v{t)+O{hn\\vm. 

Proof. Let tp{r) = 1 for \r\ < 1. Then for any e,S > 0, 

Op ^(V(p(0, X, 0/{e^h^-^')))vo = vo + 0{h^), 
since p(0,x,^) is comparable to |(a;,^)p. Let I{t) be the forward propagator for P: 

({hDt + P{t,x,hD))I{t) =0, 

jj(0) = id L2^L-- 

Let ^' = Oph{ip{p{0,x,0/{<^^h^'^^))), and set 

m = imnty'. 

We have hDtT — [P{t,x, hD),T{t, x, hD)], and by Egorov's theorem, WFft,r is contained in the flowout 
by exp{tHp) of {p(0,x,^) < 2€^h^~^^}. Now the flowout of exp(ti7p) no longer preserves the level set of p 
because p depends on t, however, if {x{t),^(t)) is an integral curve, then 

j^pit,xit),m)^ptit,x{t),at)), 

so that there is a constant C > 0, independent of h so that 

~Cp<^^p{t,x{t),at))<Cp, 
by the homogeneity of p for (x,^) in a neighbourhood of (0,0). Hence there is a constant cq so that 

Co-V(0,a;(0),e(0)) <p(t,.T(t),e(t)) < cop(0, x(0), ^(0)) 

on the flowout for < t < T. Hence a neighbourhood of (0, 0) of order h^^^^^ stays of the same order, 
although the size of e > may increase. 

We have yet to show the asserted identity property acting on v. But for this, we simply note that 
v{t) = I{t)vo, and I{t)* = /(i)^^ to write 

rv{t) = I{t)'^vo = I{t){v„ + 0{h'^)L2) = v{t) + I{t)0{h°°)L2, 

and since /(<) is unitary, we have proved the Lemma. 

□ 
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We now know that modes and quasimodes are concentrated on a scale of h^/"^^^ for any 5 > Q. We can 
measure the distance to T in the transversal direction at the point t on F using x, so it is convenient 
to have cutoffs which in a sense depend only on p, and more specifically nearly commute with P{t, x, hD). 

Lemma B.2. Let P{t, x, hD) be as above and fix N > 0. Fix e > sufficiently small and fix < a < b < e, 
and suppose ifP G C^(M) has support suppiys" C [a, b]. Then for each S > and for each I < j < N , there 
exist symbols ip^ G Si/2-s such that 

i If (f°{t,x,^,h) = (p"{p{t,x,^)/h^~'^^), then (fP € Si/2-s and ifj^ = Op/i((^°) satisfies 

[P,^P°] = 0{h'/-'+')L^^L^, 

h snppip^ C {a/i^/^""^ <P< n supp VVi", for each I < j < N, and 

iii if ip-' = Ophi'P''), then 

N 

satisfies 

with R : bounded with compact h-wavefront set. 

Proof. The proof is relatively standard, however the addition of the periodic "boundary conditions" in t 
adds a small difficulty, so we reproduce the basic argument here. 

The principal symbol of P is 

p{t,x,^) = a'^'(t,x)C,0 +b'^t)x,Xj, 

with both a*-', b^^ positive definite matrices. By the homogeneity of the quadratic forms, we clearly have 
G Si/2-s- The symbol calculus then implies the assertion (i), once we observe that -0" cuts off to a 
compact region, and hence p is bounded there. 

We now proceed to construct the ip^ satisfying (ii)-(iii). For notation simplicity, denote 

A = < p < bh^/^-^} n supp V^A° C T*M. 

From (i) and the symbol calculus, we have 

[P,V''] = ^Op,({p,/}) + i?o, 

where Rq ^ h^+'^^ Op ,i{ro) + Ro.n, where tq € Si/2-s with supprg C A, and Rq^n = 0(/i"^^/^+*'')l2->l2 
has compact ft,-wavefront set. We now compute for an arbitrary choice of ip^ G Si/2-5 with compact 
support: 

where i?j ~ h^^^^ Op h{rj) + Rj^N- As in the case of (^'^, here G 5i/2-5 has support contained in supp ViyS-' 
and Rj,N — Oi2_^^2 has compact /i-wavefront set. 

We observe that {p, ip^} has a prcfactor of since tp^ G iSi/2-5 but p G So, so we want to construct 

p^ so that iy9^}/z cancels the term of order Assume for 1 < fc < j — 1 we have 

found p^ satisfying (ii) such that 

with the rfe_i and the error in 5i/2-5 with support in A. Then if F^^^ = Sfe=o /I'^'-^/^^'^-'i/''^ satisfies 
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where R^^i = h^+^^Op hifj~i) + R]~i,N with rj_i G Si/2-s, supprj_i C A, and Rj-i,N ^ 0{h'^^'^/^+^'>) 
with compact ft,-wavefront set. Then we want to solve 

If r-,-1.0 is the principal symbol of we apply the Frobenius theorem to find such a ip^ . The support 
properties follow from the assumed support properties on rj-i and the observation that we can always 
multiply ip^ by a function of p to ensure it is supported in A, at the expense of another compactly supported 
error of order h^/^'^^. The symbolic properties follow from the assumed symbolic properties of rj^i. 

□ 

Appendix C. Real principal type 

Let us quickly show that our local principal symbol near F can be glued into a symbol of real principal 
type. Let us recall that a symbol p is of real principal type if the principal symbol is real valued, smooth, 
has compact level sets, is elliptic outside a compact set, and dp 7^ on {p = 0}. Our local model near F is 
of the form 

Pi{t,x,^) = -T + a'^{t,x)ii£,j + b'^t)xiXj, 
a*^ and b^^ are positive definite symmetric matrices, t E §>^ and r is the dual variable to t. Let pq = 
a^^ {t,x)£,i^j + b^^{t)xiXj, which is a suitable measure of the distance squared to F. Fix 5 > 0, and let 
xir) G C^(R) be equal to 1 for \r\ < S, and x(r) = for \r\ > 26. Let q ^ +P0; which is elliptic outside 
a compact set, and set 

P2 x{Po)pi + (1 - x{po))q- 
The function p2 satisfies all the requirements of real principal type, once we show that dp2 7^ on {p2 = 0}. 
For this, first note that 

{P2 - 0} C ({r = po} n {po < 26}) U {{q = 0} n {po > 6}). 

The latter set is empty, but the first is not. We observe that 

dp2 = (1 + x'iPoK'T - T^))dpo + (2(1 - x)r - x)dT, 

and that dpo ~ only if a; = ^ = 0. Now % is a non-increasing function of pq, so on {t = pq} D {po < 26}, 
x'(po)(— T — T^) > 0. Hence if {x, ^) ^ (0, 0), dp2 7^ 0. If x = ^ = 0, we have po = 0, so on {t = po} n {po < 
2(5}, T = also. But 

{2{l~x)r~x)dr^0 
for r in a neighbourhood of 0, which shows dp2 7^ on {p2 = 0}. 

Appendix D. Rescaled wavefront sets: an example 

In this section, we provide for the reader's convenience an example where one encounters the rescaled 
wavefront sets. Let us consider the quantum harmonic oscillator 

P^{hD,y+x\ 

The eigenfunctions Puj = EjUj, are semiclassical Hermite polynomials times semiclassical Gaussians. If 
uo{x) = co/i~^/'^e~^ /^'^ is the (L^-normalized) semiclassical Gaussian, and if xi^) C^(l^) is equal to 1 
near 0, then clearly 

\\{l-x{h-'x))uo\\^0{h°^) 

in any seminorm, provided 6 < 1/2. On the other hand, a simple computation shows that the semiclassical 
Fourier transform of uq 
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SO that Th{uo) has the same locaHzation in ^ as uq did in x. In other words, for any < (5 < 1/2, 

WF;,,i-,^K) -{(0,0)}. 
A similar argument imphes the same is true for any Uj provided Ej = 0{h). 

On the other hand, let us see how to use Corollary 14. 41 to prove the same localization. If Ej = 0{h), 
then the principal of P — Ej is p = + . By homogeneity, we can rescale 

which is a symbol in Sg'~^^ . As this symbol is 0{h°°) only at (x,^) = (0,0), we get the same result as 
above. 

Appendix E. Harmonic oscillator eigenfunctions 

Let Pq = —(P/dx'^ + x'^ be the one-dimensional quantum harmonic oscillator operator. The theory of 
the eigenfunctions of Pq is well established, however we need several important facts recalled here. 

Proposition E.l. Each eigenfunction it„, n = 0, 1,2, . . of Pq is a polynomial times a Gaussian: 

u„{x) = i^„(a;)e-^'/^ 

with eigenvalue A„ = 2n + 1. 

The Un can be taken to he real-valued, they form a complete orthonormal basis of (R) , and the zeros 
of Un are simple. 

Proof. The proof as usual is by creation and annihiliation operators. Let 

A± ^ Dx± ix, 

so that A± = A'^ and 

A+A- = Po - 1 = A-A+ - 2. 

A computation shows 

A_e-^'/2=0, 

so that 

A+A_e-^'/2 = = (Po - l)e-"'/2. 
The other eigenfunctions are constructed with the creation operator: 

n — 4"p-^^/2 

and a simple computation shows the u„ is a polynomial times the Gaussian e~^^/^. The creation operator 
A^ = —i{dx — x) is —i times a real-valued operator, so the polynomials can be taken to be real-valued. By 
induction, we have 

Poun = {A+A^ + l)u„ 

= A+{A^Al)uo + Un 
= A+{Pa + l)Al-\o + Un 
= A+{2{n-l) + 2)u„-i+Un 
= {2n + 1)m„. 

A simple computation shows the eigenfunctions form a complete orthnormal set. 
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To show the zeros are simple, we again assume for induction that u„ has only simple zeros, and suppose 
Un+i{xo) = 0. Then 

A^U„+l = A^A+Un 

= (Po + l)Wn 

= {2n + 2)u„, 

so 

"^'n+li^) + XUn+l{x) = i{2n + 2)Un{x), 

and if Un+i{xo) = u'j^^i(xq) — 0, then xq is a zero of u„ as well. Differentiating again, and using the fact 
that Un+i is an eigenfunction, we get 

i{2n + 2)u'^{x) = +u„+i(a:;) +xu'„+i(a;) 

= -(2?! + 2)u„+i(.t) +X^Un+l{x) + xu'^^-^^{x), 

so if Un+i{xo) ~ u'j^_^_i{xo) = 0, then u'^{xo) = as well, which contradicts the induction hypothesis. Hence 
Un+i{x) has only simple zeros. 

□ 

We are interested in these properties of the quantum harmonic oscillator eigenfunctions because, for the 
nonlinear problem studied in this note, we will take a non-smooth function of these eigenfunctions, and we 
want to understand the singularities. Let p > 0, let Un{x) be an eigenfunction of the quantum harmonic 
oscillator, and set 

V{x) = \Un{x)\^UrXx). 

Proposition E.2. The function v(x) is rapidly decaying, v(x) G H , and 

where Sci '^~^' is the space of classical symbols of order —2 — p. 

In particular, we are interested in semiclassical rescaling, and to what extent v is localized in phase 
space. The function v is not smooth, so it does not have compact semiclassical wavefront set, but because 
of the symbolic assertion in the previous proposition, there is some decay at infinity, as described in the 
next corollary. 

Corollary E.3. For any (5 > 0, the function v{x) satisfies: 

Hx)\\HiiBi0M-^)<=)=O{hn. 

Moreover, for any < 7 < 1 and < s < 3/2, the semiclassical Fourier transform satisfies 

In particular, if xi^iC) G C^(R^) is 1 in a neighbourhood o/(0, 0) , then 

x{h^x,h^-''D^)v = v + E, 

where for any < s < 3/2, 
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Proof of Proposition \E.S[ Each zero of u„ is a simple zero, so we can write 

n 
1=0 

with vq €z S, and for 1 < I < n, vi has compact support containing a single zero of «„• If xi is a zero of ii„ 
contained in the support of vi , then 

vi{x + xi) ^ \xY'x 

near x = 0, and vi{x + xi) is smooth and compactly supported away from from a; = 0. Then vi{x + xi) is 
conormal at x = 0, which implies vi is a symbol in class 5ci^ Summing in /, and using that I'o G 5ci °°, 
we get V G 5ci^^ as claimed. 

□ 
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